Ab initio two-dimensional multiband low-energy models of EtMe3Sb[Pd(dmit)2]2 and 
k-(BEDT-TTF) 2 Cu(NCS) 2 with comparisons to single-band models 
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We present ab initio two-dimensional extended Hubbard-type multiband models for 
EtMe 3 Sb[Pd(dmit) 2 ]2 and k-(BEDT-TTF) 2 Cu(NCS) 2 , after a downfolding scheme based on the 
constrained random phase approximation (cRPA) and maximally-localized Wannier orbitals, to- 
gether with the dimensional downfolding. In the Pd(dmit)2 salt, the antibonding state of the 
highest occupied molecular orbital (HOMO) and the bonding/antibonding states of the lowest un- 
occupied molecular orbital (LUMO) are considered as the orbital degrees of freedom, while, in the 
k-BEDT-TTF salt, the HOMO-antibonding/bonding states are considered. Accordingly, a three- 
band model for the Pd(dmit) 2 salt and a two-band model for the k-(BEDT-TTF) salt are derived. 
We derive single band models for the HOMO-antibonding state for both of the compounds as well. 
The HOMO antibonding band of the Pd(dmit) 2 salt has a triangular structure of the transfers with 
a one-dimensional anisotropy in contrast to the nearly equilateral triangular structure predicted in 
the extended Hiickel results. The ratio of the larger interchain transfer tj, to the intrachain trans- 
fer t a is around ti/t a ~ 0.82. Our calculated screened onsite interaction U and the largest offsite 
interaction V are ~0.7 eV and ~0.23 eV, respectively, for EtMe3Sb[Pd(dmit) 2 ] 2 and ~0.8 eV and 
~0.2 eV for k-(BEDT-TTF) 2 Cu(NCS) 2 . These values are large enough compared to transfers t 
as ~55 meV for the Pd(dmit) 2 salt and ~65 meV for the k-BEDT-TTF one, and the resulting 
large correlation strength (U—V)/t ~ 10 indicates that the present compounds are classified as 
the strongly correlated electron systems. In addition, the validity whether the present multiband 
model can be reduced to the single-band model for the HOMO-antibonding state, widely accepted 
in the literature, is discussed. For this purpose, we estimated the order of vertex corrections ignored 
in the cRPA downfolding to the single band model, which is given by W'/D, where W' is a full- 
screened-interaction matrix element between the HOMO-antibonding and other bands away from 
the fermi level (namely HOMO-bonding or LUMO-bonding/antibonding bands), whereas D is the 
energy distance between the fermi level and the bands away from the fermi level. In the present 
materials, W'/D estimated as 0.3-0.5 signals a substantial correction and thus the exchange process 
between the low-energy HOMO-antibonding and other bands away from the fermi level may play 
a key role to the low-energy ground state. This supports that the minimal models to describe the 
low-energy phenomena of the organic compounds are the multiband models and may not be reduced 
to the single-band model. 

PACS numbers: 71. 15. Mb, 71.10.Fd, 71.20.Rv, 74.70.Kn 



I. INTRODUCTION 

Searching for possible candidates of a quantum spin 
liquid state has been one of the central issues in 
condensed matter physics. Organic materials provide 
an important research areai^ and a number of or- 
ganic insulators, especially k-(BEDT-TTF) 2 Cu 2 (CN) 3 
(Ref. 0) [where BEDT-TTF is bis(ethylenedithioV 
tetrathiafulvalene] and EtMe 3 Sb[Pd(dmit) 2 ]2 (Ref. |j| 
(where dmit is l,3-dithiole-2-thione-4,5-dithiolate) are 
strong candidates of the realizations of the spin liquids. 
Although the nearest neighbor antifcrromagnctic cou- 
plings arc J^220-250 K for both compounds as specu- 
lated from the high-temperature magnetic susceptibility 
dataj 3 -^ they do not exhibit long-range magnetic order 
at least down to the temperature T^5 K. The mech- 
anism of survival of the spin degree of freedom at low 
temperatures is yet to be clarified. Strong electronic cor- 
relation and the resulting large quantum fluctuation are 



proposed to realize such a ground stated— especially for 
low-dimensional systems, but it is not so simple to estab- 
lish and identify the real materials as spin liquids. The 
low-energy structure of these organic compounds is often 
described in terms of a single-band model crossing the 
fermi level. This corresponds to the dimerization limit of 
the anti-bonding band of the highest occupied molecular 
orbital (HOMO) of organic molecules forming the dimer. 
In this limit, the system is described by a triangular- 
lattice structure and the geometrical frustration in the 
lattice is considered as a possible origin of the quantum 
fluctuation. On the other hand, however, an apprecia- 
ble anisotropy of the transfer structure in the triangu- 
lar lattice has been reported by several ab initio density 
functional calculations^— In addition, a recent dielec- 
tric measurement^ suggests that the system exhibits a 
relaxer-type response at low temperatures, T^60 K, thus 
indicating that the charge degree of freedom might bring 
about some kinds of disorder in the low-energy electronic 
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state. This is clearly beyond single-band physics. 

Recently, the ground-state properties of 
EtMe 3 Sb[Pd(dmit) 2 ]2 arc of great interest, for which 
the extremely low-tempcraturc measurements have been 
performed. 13 C NMR spectra under 7.65 T do not show 
significant broadening down to 19.4 mK^ while the 
spin- lattice relaxation rate T^ 1 exhibits a sharp drop be- 
low 1 K with a temperature dependence cc T 2 . It might 
suggest a continuous phase transition around IK due to 
some unknown symmetry breaking accompanied by the 
existence of gapless excitations in the low-temperature 
phase. On the other hand, the temperature dependence 
of heat-capacity-^ and thermal-conductivity^ data give 
a T-linear term down to O.IK, in contradiction with 
the T 2 behavior expected from . Furthermore, 
the magnetic torque measurement indicates a Pauli 
paramagnetic-like uniform magnetic susceptibility down 
to 30 niKjiS indicating a feature similar to the one- 
dimensional Mott insulator represented by Heisenberg 
or Hubbard models at half filling at low temperatures. 
Such seemingly contradicting experimental data require 
further careful theoretical analyses on the origin and 
mechanism of the spin liquid behavior. 

The low-energy electronic structure of 
EtMc3Sb[Pd(dmit)2]2 is known to be rather unique. In 
this material, the level difference between HOMO and 
lowest unoccupied molecular orbital (LUMO) is small 
compared to the dimerization gap corresponding to the 
bonding-antibonding splitting of each orbitali^r— and 
the resulting band structure shows unusual level inver- 
sion of a pair of the HOMO-LUMO bands^^ Then 
from the lower energy, the HOMO-bonding (HOMO- 
b), LUMO-bonding (LUMO-b), HOMO-antibonding 
(HOMO-ab), and LUMO-antibonding (LUMO-ab) are 
stacked separated by a finite gap. The HOMO-ab band 
is half filled and is sandwiched by the LUMO-ab and 
LUMO-b bands residing above and below 0.5 eV of the 
fermi level, respectively [see Figflja)]. Such a unique 
band structure may easily provide potential relevance 
of multiband physics, for instance, charge fluctuations 
within a dimer corresponding to the polarization between 
the crossed HOMO-LUMO level. In addition, a valence- 
bond-solid phase exists in similar compounds*^— 
where superconductivity appears in the vicinity. On the 
other hand, the geometrical frustration in the half-filled 
HOMO-ab band)ii2iii2I as well as interactions at a dimer 
site and/or between dimer sites, can also be relevant to 
the survival of the spin-liquid phase. To clarify roles 
of various factors in the low-energy physics, ab initio 
derivation of the low-energy model from first principles 
is highly desired. 

For ab initio derivations for the low-energy model, 
a scheme based on the constrained random-phase ap- 
proximation (cRPA) (RefslH and H^) for constructed 
maximally localized Wannier orbitals (MLWO) fRefl30h 
has been applied to a wide range of materials ^ Fur- 
thermore, by utilizing a quasi-low-dimensional charac- 
ter of materials as organic conductors, the method- 



ology for obtaining ab initio effective models in re- 
duced dimensions has been developed^ and applied to 
a two-dimensional (2D) effective-model derivation for 
k-(BEDT-TTF) 2 Cu(NCS) 2 ^ In Ref. a single-band 
model for the HOMO-ab band crossing the fermi level 
was derived and analyzed by a multi- variable variational 
Monte-Carlo method to demonstrate a quantitative re- 
liability of the derived model parameters through the 
study for the metal-insulator transition. It was found 
that, while this material is experimentally located in the 
metallic region with the superconductivity at low tem- 
peratures close to the border of the metal-insulator tran- 
sition^^ the derived ab initio model and its Monte- 
Carlo solution predicted an antifcrromagnctic insulator. 
With the 20% reduction of the interaction parameters, 
the ab initio model just gave a metallic solution. As a 
possible origin of the discrepancy between the theory and 
experiment, the multiband nature missing in the model 
construction was discussed, especially in terms of dynam- 
ical effects between the target band and the closest high- 
energy band (in this material, the HOMO-b band). In 
fact, in the study on the onsite Hubbard model for the k- 
BEDT-TTF system^ with the fluctuation exchange ap- 
proximation, the symmetry of the superconducting gap 
function seems to be sensitive depending on whether the 
HOMO-b band is included in the degree of the freedom 
of the effective model or not. 

In the present paper, we derive 2D multiband models 
of EtMe 3 Sb[Pd(dmit) 2 ]2 and k-(BEDT-TTF) 2 Cu(NCS) 2 
from first principles. The single-band models are also 
derived and compared with the multiband models to un- 
derstand the essence of the multiband nature. We found 
that, for EtMe3Sb[Pd(dmit)2]2, onsite intraorbital U and 
interorbital U' interactions are both ~0.7 eV and the 
Hund's rule coupling J is 0.1-0.2 eV. An offsitc interac- 
tion V is ~0.25 eV and a transfer t is 40-50 meV. As a 
result, our estimated correlation strength {U—V)/t gives 
a substantially large value as ~10. The other compound 
k-(BEDT-TTF) 2 Cu(NCS) 2 also exhibits large correla- 
tion strength as ~10. We also estimate the reliability 
of the downfolding whether the multiband model can be 
reduced to the single-band model with the ab initio cRPA 
framework. The measure of the reliability is given by an 
order estimate of the vertex correction to the effective in- 
teraction. Wc found that, for both EtMe 3 Sb[Pd(dmit)2]2 
and k-(BEDT-TTF) 2 Cu(NCS) 2 , the vertex correction is 
not small and may well affect the low-energy physics, 
thus suggesting that the multiband analysis beyond the 
single-band one is needed for describing the low-energy 
physics properly. 

The present paper is organized as follows: In Sec. 
II, we define effective models to be derived and de- 
scribe a scheme for the derivation. Computational de- 
tails and results for EtMe 3 Sb[Pd(dmit) 2 ]2 and k-(BEDT- 
TTF) 2 Cu(NCS) 2 are given in Sec. III. We discuss charac- 
teristic aspects of the derived effective models including 
correlation strength and the estimate of the reliability of 
the downfolding in Sec. IV. Summary is given in Sec.V. 
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II. MODEL 

Here, we describe basic procedures for the derivation 
of effective low-energy models for the two compounds. 
The basis of the Hamiltonian is the Wannier function 
associated with anti-bonding/bonding states of HOMO 
or LUMO of dmit or BEDT-TTF molecules that form a 
dimer. The derived Hamiltonian is explicitly given in a 
form of the 2D multiband extended Hubbard model as 

^ = ^Z e iH a ifj,a a if*c + ^iM^ a \p,a a m^ 

(rip, a i^j p 



9 ^Mf a i^a a iup a 3vp a ina 



ap ij pv 



ap ij pv 



with (flipa) being a creation (annihilation) operator 
of an electron with spin a in the /ith Wannier orbital 
localized at the ith dimer site. The tip and tipjp param- 
eters are given by 



and 



tip = (<t>iij\H--Ks\<t>in) 



(0j//|'Hks|0j>), 
,t 



(2) 



(3) 



respectively, with \<t>ip)=a ill |0) and 'Hks being an ef- 
fective one-body Kohn-Sham Hamiltonian. In this 
model for EtMe 3 Sb[Pd(dmit) 2 ]2, we take three Wan- 
nier orbitals, i.e., those of LUMO-b, HOMO-ab, and 
LUMO-ab, to represent the Hamiltonian. In k-(BEDT- 
TTF) 2 Cu(NCS) 2 , we take two Wannier orbitals; HOMO- 
ab and HOMO-b. Note that the Wannier orbitals are 
made for each band individually, since the bands sepa- 
rate with each other. With this construction of the Wan- 
nier function, the transfers for different orbitals are zero 
(i.e., ti^j V =Q for \x^v). The Vipj V and Jipj V parameters 
in Eq. ([1} are screened Coulomb and exchange integrals 
in the Wannier-orbital basis, respectively, expressed as 

^=//^r^(r)0 4 ^r)^(r,rO^(r')0^(r')(4) 
and 

J iH»=JJ rfrdr^(r)^(r)W/(r,r')^(r')0 v (r')(5) 

with W(r, r') being a 2D screened Coulomb interaction 
in the low-frequency limit. Here, is nothing but the 
onsite Hubbard U, while Jipip is set to zero by defini- 
tion. If one keeps only the HOMO-ab band crossing the 
fermi level and treats other bands as screening bands, the 
multi-band model is reduced to a single-band model 

a ijkj ap ij 

+ 2 X) X! Jl 3 ( a L a ]p a iP a 3° + a \a a \p a iP a 3o) > ( 6 ) 



ap i^ij 



where we drop the orbital index ji. 

The calculation for W(r,r') in Eqs. (g]) and © fol- 
lows Ref. H3, where a new cRPA framework was devel- 
oped for the purpose to derive effective interactions of 
models defined in lower spatial dimensions. This new 
scheme is suitable for quasi-low-dimensional materials 
as the present system. The cRPA method is originally 
formulated in the RPA framework with the constraint 
for the band degree of freedom to eliminate only the 
degrees of freedom far from the fermi level in energy. 
This is called the band downfolding. In the proposed 
scheme of the supplementary downfolding, however, the 
concept of the constraint is additionally relaxed to in- 
clude the screening by the polarization in the other lay- 
ers/chains even within the target bands. This utilizes 
the real-space representation of the polarization function 
and elimination of specific-polarization blocks associated 
with the spatial dimensions perpendicular to the target 
layer/chain, which leads to a low-dimensional effective in- 
teraction for the target layer/chain. We call it the dimen- 
sional downfolding. Practically, the band+dimensional 
downfolding is performed in two steps: We first per- 
form the band downfolding to derive a three-dimensional 
model for a small number of bands near the fermi levels 
This is followed by the dimensional downfolding in the 
second step^ With this idea, we can naturally derive 
the low-energy model in any dimension. In the present 
case, we use it for the derivation of multiband and 
single-band 2D models for EtMc3Sb[Pd(dmit)2]2 and k- 
(BEDT-TTF) 2 Cu(NCS) 2 . 



III. RESULTS 

A. Computational Detail 

Ab initio density- functional calculations were per- 
formed by the program package xTAPP, which is a mas- 
sively parallelized version of Tokyo Ab initio Program 
Package (TAPP) (Ref. HI) and is able to perform the 
present large scale calculations efficiently. The xTAPP 
program adopts plane-wave basis sets and the present 
calculations were performed with norm-conserving pseu- 
dopotentials^ and generalized gradient approximation 
(GGA) (Ref. l39h for the exchange-correlation potential. 
The experimental structure of EtMe3Sb[Pd(dmit) 2 ] 2 at 
4 K was taken from x-ray crystallography^ while the 
atomic structure of k-(BEDT-TTF) 2 Cu(NCS) 2 at 15 K 
was taken from neutron data42 For both of the sys- 
tems, the positions of the hydrogen atoms were re- 
laxed. The cutoff energies in wavefunctions and charge 
densities were set to 36 Ry and 256 Ry, respectively. 
5x5x3 and 5x5x5 fc-point samplings were employed for 
EtMe 3 Sb[Pd(dmit) 2 ] 2 and k-(BEDT-TTF) 2 Cu(NCS]2, 
respectively. The construction of MLWO follows Ref. |3fj. 
The polarization function was expanded in plane waves 
with an energy cutoff of 5 Ry and the total number 
of bands considered in the polarization calculation was 



set to 2000 for EtMe 3 Sb[Pd(dmit) 2 ]2 and 750 for k- 
(BEDT-TTF) 2 Cu(NCS) 2 . This condition corresponds 
to considering the excitation up to ~20 eV above the 
fermi level. The Brillouin-zone (BZ) integral on the 
wavevector in the polarization calculation was evaluated 
by the generalized tetrahedron method.— In the dimen- 
sional downfolding procedure, the calculated 2D effective 
Coulomb/exchange integrals were extrapolated to values 
with the infinite number of screening layersi 33 ' 34 



B. EtMe 3 Sb[Pd(dmit) 2 ]2 

Figure [TJ (a) shows our calculated GGA band struc- 
tures (solid lines) of EtMe 3 Sb[Pd(dmit) 2 ] 2 . The band in 
[-0.1 eV: 0.3 eV] is the HOMO-ab band which is sand- 
wiched by the LUMO-b and LUMO-ab bands at -±0.5 
eV. The dotted lines represent the tight-binding band ob- 
tained with the transfer parameters in Table HI where the 
definition of each transfer is given in the panel (b). The 
transfer parameters were estimated from the MLWOs. In 
the panels (c) and (d) , we display our calculated MLWOs 
for the HOMO-ab and LUMO-b bands, respectively. As 
the initial guess for the HOMO-ab Wannier orbital, we 
used the p-type Gaussian with the width 1.67 A, cen- 
tered at the dimer center. Also, the initial guess for the 
LUMO-b and LUMO-ab MLWOs are superposition of the 
p-type Gaussians with width 1.18 A put on the S atoms 
adjacent to Pd in each dmit molecule. 



TABLE I: List of transfer parameters in the extended Hub- 
bard model of EtMe3Sb[Pd(dmit)2]2- For the definition of 
the transfers, see Fig. QJb) and e is the single-particle level 
energy. The unit is meV. 





e 


t a 


t b 


t c 


LUMO-b 


-490 





8.4 


-16.1 


HOMO-ab 


34.8 


54.4 


44.9 


40.2 


LUMO-ab 


495 


-25.7 


24.8 


15.1 


HOMO-ab exHiickel 




28.2 


27.0 


25.0 



The transfer parameters for the HOMO-ab band in 
the extended Hiickcl method are given in the bottom of 
the tabled Note that, the notation for the transfers in 
Ref. Hi] is different from the present paper; in the corre- 
spondence between the former and latter, ts=t a t a —tb, 
and t r =t c . In comparison with the ab initio values, the 
Hiickel values are smaller than the ab initio ones by the 
factor of 1/2-2/3. On top of that, the ab initio param- 
eters exhibit an appreciable one-dimensional anisotropy 
along the t a direction, i.e., t a >tb~t c . In contrast, the 
Hiickel parameters is close to the equilateral triangular 
lattice as t a ~tb~t c . 

We next show in Fig. [2] our calculated interaction pa- 
rameters for the 2D single-band model. The panels (a) 
and (b) give the diagonal and off-diagonal parts of the 
interaction, respectively. In the presentation, we define a 
reference site as a site indicated by the arrow. The value 
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(c) 
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FIG. 1: (Color online) (a) Our calculated GGA band 
structure (solid lines) and the Wannier interpolated bands 
(dotted lines) with the transfers listed in Table [T] of 
EtMe3Sb[Pd(dmit)2]2- The crystal structure consists of al- 
ternating layers (parallel to the ab plane) of the Pd(dmit)2 
anion and polymeric EtMe3Sb + cation. Band dispersions are 
plotted along the high-symmetry points in the ab plane, where 
T=(0, 0, 0), X=(a*/2, 0, 0), M=(a*/2, b* /2, 0), Y=(0, b*/2, 
0). The zero of energy is the fermi level. The HOMO-ab 
band is drawn by the dark blue dotted lines and the LUMO- 
b and LUMO-ab bands are described by light green dotted 
lines, (b) Triangular lattice and the definition of transfers. 
The system contains two equivalent Pd(dmit)2 layers and the 
dimer units stack along the a+b direction in one layer and 
the a— b direction in the other layer. In this figure, the latter- 
layer configuration is depicted. The t a direction is defined as 
the stacking direction and the tj, direction is defined to be 
parallel to the b axis. The t c direction is given as the direc- 
tion of the remaining side of the triangular lattice, (c) Our 
calculated Wannier function of the HOMO-ab state in view 
along the b axis, drawn by VESTAL The dark blue surfaces 
indicate positive isosurface and the light yellow surfaces indi- 
cate negative isosurface. (d) The same plot for the LUMO-b 
state. 



in the ellipsoid represents the interaction value between 
the reference site and this ellipsoid one, in the unit of 
eV. For example, in the panel (a), 0.61 (eV) in the ref- 
erence site is the onsite Hubbard U and 0.22 (eV) next 
on the right is the offsite Coulomb integral between the 
reference site and this site. In the figure, we show the in- 
teraction value up to the tenth neighbors and the further 
distant interactions arc negligibly small. Note that the 
system has a C 2 symmetry. The offsite exchange integral 
in the panel (b) is given in the unit of meV. In princi- 
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pie, the interaction range is short as up to the nearest 
shell and the values are so small as nearly half of those 
of k-(BEDT-TTF) 2 Cu(NCS) 2 .^ 

(a) Coulomb (eV) 
(MM IKI6 — (0.05 



(a) Coulomb (eV) 

0.06- (ins 
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\ 1/ 
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0.27, 





\ 



'0.22 


0.22 


0.23 A 
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0.23 
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0.23, 
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0.20 0.21 




0.06J 




0.68 


0.71 A 




'0.27 


0.27 


0.27" 


0.69 


0.69 




0.27 


0.27 


0.27 


0.69 


0.72 j 




v 0.27 


0.27 


0.27 



FIG. 2: Effective Coulomb (a) and exchange (b) integrals of 
the single band model for EtMe 3 Sb[Pd(dmit) 2 ]2 in Eq. (©. 
The arrow indicates a reference site. The view of the trian- 
gular lattice follows Fig. QJb) . 



The interaction values of the 2D three-band model are 
displayed in Fig. [3J In the multi-band case, the interac- 
tion is given as the 3x3 matrix as 



Voijl Vo\j2 K)lj3\ / Jbljl Joijl Joij3 

Vo2jl Vo2j2 Vo2j3 and I Jo2jl Jo2j2 Jo2j3 
V()3jl Vozj2 Vo3j3/ \Jo3jl Jo3j2 Jo3j3> 



, (7) 



where Vb w v {Jo^ju) indicates the Coulomb (exchange) 
integral between the /xth orbital in the 0th reference site 
and the isth orbital in the jth. site. In the orbital labels, 
(1,2,3) := (LUMO-b, HOMO-ab, LUMO-ab). The on- 
site Coulomb interactions are nearly ^0.7 eV and hold 
U~U'. On the offsite interactions, the orbital depen- 
dence is negligible, thus being well described as the func- 
tion of the distance between the centers of the Wannier 
functions. On the basis of this aspect, for the interactions 
beyond the first shell, we show the averaged value over 
the orbitals. The Hund's rule coupling at the reference 
site is rather large (e.g., the exchange integral between 
the LUMO-b and LUMO-ab orbitals as ~0.22 eV). The 
exchange integrals with the neighboring site are rather 
orbital dependent, and there exist rather large values 
more than 10 meV. 



(b) Exchange (meV) 



fl.22 1.78 
1.78 1.90 
1.10 1.42 
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FIG. 3: Effective Coulomb (a) and exchange (b) integrals 
of the three-band model for EtMe 3 Sb[Pd(dmit) 2 ]2 in Eq. QJ. 
The arrow indicates a reference site. The view of the trianglar 
lattice follows Fig. []Jb). 



(d), respectively The extended Hiickel transfers for the 
HOMO-ab band 3 -^ is given in the bottom of TableHU In 
k-(BEDT-TTF) 2 Cu(NCS) 2 , t a ~t b holds and the lattice 
is approximated as an isosceles triangle with t=(t a +tb)/2 
and t'=t c - For both the ab initio and extended Hiickel 
results, this trend well holds, but the former frustration 
strength t'/t as ^0.66 is somewhat smaller than the lat- 
ter one as ^0.84. We also note that, in the ab initio 
results, the tj transfer has a finite magnitude. 



C. k-(BEDT-TTF) 2 Cu(NCS) 2 

We next move to k-(BEDT-TTF) 2 Cu(NCS) 2 . Fig- 
ure S] (a) displays our calculated GGA band (solid lines) 
and the tight-binding band (dotted lines) with the trans- 
fer integrals in Table ITT] The definition of each transfer 
follows the panel (b) and the Wannier function of the 
HOMO-b and HOMO-ab are plotted in the panels (c) and 



TABLE II: List of transfer parameters in the extended Hub- 
bard model of k-(BEDT-TTF) 2 Cu(NCS) 2 . For the definition 
of the transfers, see Fig. f3Jb). The unit is meV. 





e 


ta 


t b 


t c 


t d 


HOMO-ab 


0.02 


-64.8 


-69.3 


44.2 


-11.5 


HOMO-b 


-447 


-40.3 
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(c) (d) 




FIG. 4: (Color online) (a) Our calculated GGA band struc- 
ture (solid lines) and the Wannier interpolated bands (dot- 
ted lines) with the transfers listed in Table HIl of k-(BEDT- 
TTF)2Cu(NCS)2- The crystal structure contains alternat- 
ing layers (parallel to the be plane) of BEDT-TTF donor 
molecules and polymeric Cu(NCS)2~ anions. Band disper- 
sions are plotted along the high-symmetry points in the be 
plane, where F=(0, 0, 0), Y=(0, b*/2, 0), Z=(0, 0, c*/2), 
and M=(0, 6*/2, c*/2). Note that the a axis is interlayer 
axis. The zero of energy is the fermi level. The HOMO-ab 
band is drawn by the dark blue dotted lines and the HOMO-b 
band is described by light green dotted lines, (b) Definition 
of the transfers. The directions of the b and c axes forming 
the plane are also drawn, (c) Our calculated Wannier func- 
tion of the HOMO-b state in view along the a axis, drawn by 
VESTAL The dark blue surfaces indicate positive isosurface 
and the light yellow surfaces indicate negative isosurface. (d) 
The same plot for the HOMO-ab state. 



The interaction parameter for the 2D single-band 
model is given in Fig. [5j where the panels (a) and (b) 
are diagonal and off-diagonal parts for the interaction, 
respectively. How to see the figure is the same as Fig. [5] 
The interaction range for k-(BEDT-TTF) 2 Cu(NCS) 2 is 
shorter than that of EtMe3Sb[Pd(dmit)2]2 a s far as 
we compare the distance with the values up to ^0.05 
eV. This is due to the fact that the interdimer dis- 
tance of k-(BEDT-TTF) 2 Cu(NCS) 2 is larger than that of 
EtMe 3 Sb[Pd(dmit) 2 ]2 [7.68 A for the k-BEDT-TTF salt 
and 6.31 A for the Pd(dmit)2 salt]. Furthermore, in gen- 
eral, the in-plane interaction range is roughly determined 
by the interlayer distance ) 33 ' 34 because, in the dimen- 
sional downfolding, the interlayer metallic screening is 



switched on from this distance [16.44 A for the k-BEDT- 
TTF salt and 18.52 A for the Pd(dmit) 2 salt] and, within 
this range, the in-plane interaction gives finite values. 
We also mention that the exchange value of k-(BEDT- 
TTF)2Cu(NCS)2 is nearly twice as large as that of 
EtMe3Sb[Pd(dmit)2]2, which may reflect the difference in 
the dimer orientation of k-(BEDT-TTF) 2 Cu(NCS) 2 and 
EtMe 3 Sb[Pd(dmit)2] 2 [Compare Fig. [1(b) and Fig.^b)]. 



(a) Coulomb (eV) (b) Exchange (meV) 




FIG. 5: Effective Coulomb (a) and exchange (b) integrals 
of the single band model for k-(BEDT-TTF) 2 Cu(NCS) 2 in 
Eq. J6]). The arrow indicates a reference site. The view of the 
triangular lattice follows Fig. |4fb). 



We next show in Fig. [5] our calculated 2D interaction 
parameters for the two-band model. The interactions are 
given by the matrix form as 

( JWJW^) and ( J ^ J ^A . ( 8 ) 

\VQ2j\VQ2j2 ) \J02jl •J02j2 J 

In the orbital labels, (1,2) := (HOMO-ab, HOMO-b). 
The notations are the same as Fig. [31 but note that the 
symmetry of the interaction matrix is different from that 
for EtMc3Sb[Pd(dmit)2]2- The onsite Coulomb interac- 
tion is ~0.8 eV and the Hund's rule coupling is ~0.36 eV, 
being larger than those of EtMc 3 Sb[Pd(dmit) 2 ]2- The 
offsitc Coulomb integrals are not orbital dependent, while 
the offsite exchange integrals are largely orbital depen- 
dent and the largest value is as large as ^30 meV. 



IV. DISCUSSION 

Here, we first remark that the antifcrromagnctic ex- 
change coupling J is estimated in the strong-coupling 
perturbation as J - At 2 /{U~V) = 245 K for the 
Pd(dmit) 2 salt and J - 360 K for k-BEDT-TTF salt, 
when we use t = (t a + % + t c ) /3 to compare with the ex- 
perimental data deduced from the fitting to the isotropic 
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FIG. 7: (Color online) A schematic diagram of the screen- 
ing process of the onsite interaction for the HOMO-ab 
Wannier orbital of EtMe3Sb[Pd(dmit)2]2, generated by the 
band+dimensional downfolding. 



FIG. 6: Effective Coulomb (a) and exchange (b) integrals of 
the two band model for k-(BEDT-TTF) 2 Cu(NCS) 2 in Eq. JTJ. 
The arrow indicates a reference site. The view of the trian- 
gular lattice follows Fig. Ub). 



triangular Heisenberg model in the high-temperature 
expansion^ When we consider the insufficiency of the 
single-band model and the experimental uncertainty 
coming from the lack of the data above 300 K, the re- 
sults are favorably compared with 250 K estimated for 
the both compounds from such fittings of the magnetic 
susceptibility. 

Now we discuss the screening effects due to the band 
downfolding from the multiband to single-band mod- 
els and due to the dimensional downfolding as well. 
In Fig. [71 we draw a schematic diagram showing how 
the onsite interaction value of the HOMO-ab band of 
EtMc3Sb[Pd(dmit)2]2 is screened. The bare Coulomb 
interaction is 3.55 eV and in the stage of the 3D three- 
band model with the conventional band downfolding, the 
interaction value becomes 0.89 eV. With considering the 
dimensional downfolding, the value is reduced further to 
0.69 eV by 0.2 eV. When we consider the screening due 
to the LUMO-b and LUMO-ab bands in addition, to con- 
struct the single-band model, the value becomes 0.61 eV 
(the reduction of 0.08 eV). So, the largest part of the 
screening comes from the high-energy-electron screening 
considered in the initial band-downfolding stage, while 
the screening effects due to the dimensional downfold- 
ing, as well as the reduction to a single-band model may 
not be neglected quantitatively. We note that the fully- 
screened-RPA value is 0.14 eV. 

The same diagram for k-(BEDT-TTF) 2 Cu(NCS) 2 is 
shown in Fig. [5J In this material, while the reduction of 
the value due to the dimensional downfolding is nearly 
the same as the EtMc3Sb[Pd(dmit)2]2 case, the reduc- 
tion due to the band downfolding from the two-band to 
single-band model is appreciable (the reduction of ^0.2 
eV). This is probably due to the fact that the spatial 
distribution of the eliminated HOMO-b orbital is similar 
to that of the target HOMO-ab orbital and largely over- 
laps, resulting in a more efficient screening. The bare 



and fully-screened-RPA values are 3.61 and 0.19 eV, re- 
spectively. 
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FIG. 8: (Color online) A schematic diagram of change in the 
onsite interaction for the HOMO-ab Wannier orbital of k- 
(BEDT-TTF) 2 Cu(NCS) 2 , generated by the band+dimension 
downfolding. 



We next discuss the degree of the correlation strength 
of the 2D effective models, which is measured by 
(U—V)/B, where U and V are the onsite and offsitc 
Coulomb interactions, respectively, and B is the band- 
width. Since we are interested in low energy, we fo- 
cus on the HOMO-ab band. Table Mil lists our calcu- 
lated parameters and the correlation strength. For both 
EtMe 3 Sb[Pd(dmit) 2 ]2 and k-(BEDT-TTF) 2 Cu(NCS) 2 , 
(U—V)/B is comparable to the unity and thus these sys- 
tems are classified as the strongly correlated electron sys- 
tem. Although the correlation strength of the single-band 
model is somewhat smaller than that of the multiband 
model, it is still large as close to 1. 

We now detail the validity of reducing from the multi- 
band to single-band model in the present case, by apply- 
ing a general criterion^ In cRPA for deriving the effec- 
tive interaction of the single-band model, the multiband 
is divided into the target HOMO-ab band and the rest. 
Then, with excluding a polarization formed in the target 
band, the effective interaction is calculated. Figure [9] (a) 
exhibits a schematic diagram for deriving the single-band 
effective interaction illustrated symbolically in the second 
order one, where low-energy electrons in the HOMO-ab 
band (external lines) is screened by a polarization (in- 
ternal lines) via high-energy electrons (namely, electrons 
belonging to the bands away from the fermi level) in the 
eliminated bands [i.e., the LUMO-ab/LUMO-b bands for 



8 



TABLE III: Our estimated correlation strength (U—V)/B and downfolding measure W'/D. Here, U and V are onsite and 
offsitc Coulomb interactions of the HOMO-ab Wannier orbitals, respectively, and B is the bandwidth of the HOMO-ab 
band. Also, W' is the matrix element of the full-RPA interaction between the HOMO-ab and LUMO-b Wannier orbitals 
for EtMe 3 Sb[Pd(dmit) 2 ]2 and HOMO-ab and HOMO-b Wannier orbitals for k-(BEDT-TTF) 2 Cu(NCS) 2 and D is the mini- 
mum energy distance between the fermi level and bands away from the fermi level. The unit for U, V, B, W' , and D are given 
in eV. 
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the Pd(dmit) 2 salt and the HOMO-b band for the k- 
BEDT-TTF salt]. Here, the capital characters "L" and 
"H" describe propagators of the low- and high-energy 
electrons, respectively. In the figure, an electron-hole pair 
formed in the low- and high-energy bands is displayed. 
The infinite series of such a diagram gives the effective 
interaction used in the single-band model. In this treat- 
ment, the vertex which modifies the interaction is totally 
dropped. In the panels (b) and (c), we show the correc- 
tion to the original interaction process in (a); the vertex 
in the shaded area corrects the polarization [panel (b)] or 
introduces an additional interaction process [(c)]. To the 
first approximation, the vertex is approximated by the 
single fully screened RPA interaction; the correction to 
the polarization is written as an exchange process (d) and 
the additional interaction process for the propagators is 
written as a single ladder process (e). 

Let us consider the magnitude of these corrections 
to the original process. The polarization is given by 
—iGhGl with Gl and Gh being the propagators of 
the low- and high-energy electrons, respectively, and this 
magnitude roughly scales as — (1/-D), where D is the min- 
imum energy distance of the high-energy band from the 
fermi level. We note that, in this estimate, the effect of 
the transition matrix element is dropped for the simplic- 
ity. Similarly, the correction term described in the panel 
(d) is given by G h G l G h G l W L h ~ W LH /D 2 with W LH 
being the fully-screened interaction between the low- and 
high-energy electrons^ Thus, the order of the correc- 
tion to the original polarization is scaled as Wlh/D 
From the similar discussion, the order of the correction 
described in the panel (e) is also given as Wlh/D. This 
value of Wlh/D is a measure for a stable downfolding 
from a multiband to a single-band model. When this 
value is small, the vertex correction to the cRPA interac- 
tion can be neglected. On the contrary, when the value is 
large, the cRPA downfolding breaks down and one must 
solve the multiband model directly without reducing to 
the single-band model. 

We show in Table HTT1 our estimated Wlh/D, where 
we have approximated Wlh by the full RPA form W'; 
in EtMe 3 Sb[Pd(dmit) 2 ] 2 , W is calculated as the full- 
RPA Coulomb interaction matrix element between the 
HOMO-ab and LUMO-b Wannier orbitals, and, in n- 
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LL LL L L or H L 




LL LL L L or H L 



(d) (e) 




FIG. 9: (Color online) (a) A schematic diagram for effective 
interactions of the single-band model, where we draw sym- 
bolically the second-order one. In cRPA, the interaction be- 
tween low-energy electrons in the HOMO-ab band (external 
lines) is screened by polarizations (internal lines) via high- 
energy electrons, where the capital characters "L" and "H" 
describe propagators of the low- and high-energy electrons, 
respectively. In the panel (b), we show the correction to the 
effective interaction (a) , where the polarization is modified by 
the vertex (shaded area) ignoring the RPA process, while in 
(c), the vertex introduces an additional interaction process. 
In the first-order approximation, the vertex is replaced by the 
fully screened interaction. We draw the schematic diagrams 
for those; (d) the vertex correction to the polarization and (e) 
the ladder-type interaction for the propagators. 



(BEDT-TTF) 2 Cu(NCS) 2 , W' is estimated as the matrix 
element between the HOMO-ab and HOMO-b Wannier 
orbitals. The resulting W'/D is as large as 0.3-0.5. This 
value is rather large compared to the transition-metal 
oxides. In the transition-metal oxide, d bands (i 2ff or e g 
bands) form the low-energy band and the closest high- 
energy bands consist of oxygen-p bands. In this case, the 
matrix element of the full-RPA interaction Wd p between 
^ig/s-g and Op orbitals is substantially smaller, because 
the centers of the Wannier orbitals are spatially apart 
with each other and the full-RPA interaction does not 
work for such a length scale. On the other hand, in the 
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case of the organic compounds, the Wannier orbitals of 
the low- and high-energy bands have the same centers 
and therefore the matrix element for the full-RPA inter- 
action does not disappear. This is a characteristic fea- 
ture of the organic compound compared to the transition- 
metal oxide. As a result, in the case of the organic com- 
pound, the downfolding treatment based on cRPA to re- 
duce to the single-band model would be insufficient and 
the multiband analysis is strongly recommended. 

One of the compounds exhibiting quantum spin-liquid 
behavior is k-(BEDT-TTF) 2 Cu 2 (CN) 3 . The ab initio pa- 
rameters for the 3D single-band systems have been de- 
rived in Ref§. Since the HOMO bonding orbital appears 
to be entangled with other lower-energy bands, one may 
need to treat more than two bands for the multi-band 
description, in contrast to k-(BEDT-TTF) 2 Cu(NCS) 2 , 
which is left for future studies. 



V. CONCLUSION 

To summarize, in the present paper, we derived the 
2D multiband models for EtMc 3 Sb[Pd(dmit) 2 ] 2 and n- 
(BEDT-TTF) 2 Cu(NCS) 2 , with the ab initio cRPA plus 
MLWO framework suitable for strongly correlated elec- 
tron systems. The correlation strength for the derived 
model, (U—V)/B, is estimated as a large value of ~1, 
as shown in Table IIII1 indicating the strongly corre- 
lated nature of electrons in these systems. The trans- 
fer structure of EtMc3Sb[Pd(dmit) 2 ] 2 is characterized by 
a lattice of quasi-2D scalene triangular lattice with ID 
anisotropy, whereas that of k-(BEDT-TTF) 2 Cu(NCS) 2 



approximately has a 2D square-lattice structure with 
a weaker next-nearest neighbor transfer. To discuss 
whether the multiband model can be reduced to the 
single-band model with the usual cRPA treatment, we 
estimated the order of the vertex correction to the 
cRPA interaction. We found that the vertex process be- 
tween the target HOMO-ab band and others [i.e., the 
LUMO-ab/LUMO-b bands for the Pd(dmit) 2 salt and 
the HOMO-b band for the k-BEDT-TTF salt] is not 
small. Therefore, the multiband model should directly 
be analyzed, not by reducing to the single-band model, 
to clarify the realistic correspondence with strongly cor- 
related phenomena including non-fermi liquid behavior, 
quantum spin liquid phase and/or unconventional super- 
conducting mechanism found in the real organic com- 
pounds. This remains to be explored. 
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